Interacting drift-diffusion theory for photoexcited electron-hole gratings in 

semiconductor quantum wells 
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Phase-resolved transient grating spectroscopy in semiconductor quantum wells has been shown to 
be a powerful technique for measuring the electron- hole drag resistivity pe/i, which depends on the 
Coulomb interaction between the carriers. In this paper we develop the interacting drift-diffusion 
theory, from which peh can be determined, given the measured mobility of an electron-hole grating. 
From this theory we predict a cross-over from a high-excitation-density regime, in which the mobility 
has the "normal" positive value, to a low-density regime, in which Coulomb-drag dominates and 
the mobility becomes negative. At the crossover point, the mobility of the grating vanishes. 



The phenomenon of Coulomb drag, whereby an elec- 
tronic current driven in a quasi-two dimensional electron 
gas drags carriers in an adjacent quantum well, creating 
a potential difference or a current in the latter has been 
a topic of intense interest in condensed matter physics 
for the past two decades [ll-[lo|. The effect was first pre- 
dicted, theoretically, by Pogrebinskii [H and Price [1], 
but the field exploded only in the early 1990s, follow- 
ing the first successful realization of independently con- 
tacted bilayer structures Since then the field has 



flourished, due to the insights it affords on the nature of 
Coulomb correlations, non-equilibrium fluctuations, and 
quantum coherence between spatially separated systems. 
While the original experiments were done on electronic 
bilayers Q, it was soon realized that electron- hole bi- 
layers offer even more interesting scenarios 0, H, [13, [l3 ■ 
For example, electrons and holes can condense in an exci- 
tonic superfluid, with or without a magnetic field 
[l5l |. and the occurrence of such a condensation should 
lead to striking manifestations in Coulomb drag exper- 
iments [15|. More recently, drag effects have also been 
investigated for different spin populations in a single 
quantum well (sp in Coulomb drag) Til- 21 1, in bilayer 
graphene (iil. l23l|. in topological insulators 24 1, and in 
trapped cold atoms [25| . 

The study of drag effects in electron- hole bilayers is no- 
toriously complicated by the fact that, in order to host 
carriers of opposite polarities, the two layers must be held 
at different chemical potentials, while their spatial sep- 
aration is of the order of 20-30 nm. In contrast to this, 
a non-equilibrium non-homogeneous distribution of elec- 
trons and holes can be rather easily achieved in a single 
layer with a laser of frequency larger than the band gap, 
which creates an equal number of electrons and holes, in 
addition to the carriers (say electrons) that are already 
present at equilibrium. 

In a recent series of experiments [i^, 2^ 27 [ the inter- 



ference between two laser beams coming from different 
directions and polarized in the same direction has been 
exploited to create a transient electron-hole (e-h) grating, 
i.e., a spatially periodic modulation of the electron and 
hole densities on the surface of an n-type GaAs quan- 



tum well. Electron- hole drag manifests itself in quite a 
striking and direct way in the mobility of the e-h grating 



under the action of an electric field [27|, |28| . The crucial 
observation is that the drift velocity of the e-h grating - a 
collective formation - differs, in general, from the drift ve- 
locity of the background majority carriers (electrons). If 
the background carriers were absent (e.g. in an intrinsic 
material) then the drift velocity would vanish, since the 
net electric force on a neutral object is zero. In the oppo- 
site limit, in which the e-h grating is a small perturbation 
on the background electron density, the mobility depends 
crucially on the rate of momentum transfer between elec- 
trons and holes [l^. Let us assume at first that this can 
be neglected. Then the constraint of charge neutrality 
causes the dense electrons to follow the few holes rather 
than the other way around. Under these conditions the 
grating moves with the drift velocity of the holes, and 
the mobility is positive [i^. Coulomb collisions between 
counter-flowing electrons and holes produce the equiva- 
lent of an electric field that pushes the holes in the direc- 
tion opposite to the external field. Which of the two fields 
ultimately dominates depends on the mobility of the elec- 
trons: if that mobility is sufficiently high the holes will 
inevitably reverse their direction of flow. The mobility 
of the grating will then be negative and approach the 
mobility of the electrons when the latter is large. 

The possibility of observing anomalous (i.e., negative) 
mobility for an electron-hole packet in n-type semicon- 
ductors was first pointed out by McLean and Paige [s^ 
and recently confirmed by Yang et al. in n-GaAs. 
However, these authors relied for their analysis on a sim- 
plified phenomenological model in which the e-h grating 
and the background electrons are treated as separate en- 
tities. The e-h grating, being a neutral entity, does not 
directly respond to the electric field. Hence, it is impos- 
sible, within this model, to describe the competition be- 
tween Coulomb drag and the "normal" ambipolar mobil- 
ity of the grating: the resulting formulas are correct only 
in the limit of strong Colulomb drag and cannot be used 
to predict the cross-over between the "normal" and the 
"anomalous" (i.e. Coulomb drag-dominated) regimes. 

In this paper we develop a full-fledged interacting drift- 
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diffusion theory from which we derive general formulas 
for the mobility and the diffusion constant of the e-h grat- 
ing. These formulas not only enable us to extract pre- 
cise values of the Coulomb drag resistivity from Doppler 
velocimetry data, but they also lead to a striking predic- 
tion: the e-h grating mobility can be driven through a 
sign reversal by changes in excitation intensity, tempera- 
ture, or background density. With an appropriate choice 
of parameters, the e-h grating can be made stationary in 
the presence of an electric field. 

Our starting point is the drift-diffusion equation for 
electron and hole densities {n and p respectively), in 
which, however, we take into account the presence of the 
off-diagonal homogeneous conductivity aeh (this is what 
makes our theory interacting, whereas, in the "standard" 
theory, a^h is set to zero). Thus, we have 



ductance matrix can be expressed as [IE 



edt 



P 



= -V • 



(Thh 



(1) 



where the effective electric fields are expressed as E°^^ — 
E ± ^Vfj,c/y Here, E includes both external and built- 
in electric fields, /ic and fiy are chemical potentials of 
electrons and holes, respectively, not to be confused with 
the homogeneous mobilities /ig and fih, which will appear 
later. By substituting the electric fields, the continuity 
equations become 

edtn = V • [eZ^eeVn + eDeh^p] + E ■ Vag 

+ CTeV • E , (2) 

edtp - V • [eDhhyp + eDhe^n] - E • Vcr,, 

- fT,,V • E , (3) 

where Ue = cTee + Ce/i and Uh — crhh + <^eh ■ The diffusion 
matrix elements are Daj3 = {—1)^""' aaj{x~^)'yi3 /^'^ (sum 
over 7 implied), where x-yP = dn^/d^p is the suscepti- 
bility matrix 31|. 

Following the standard procedure (29l | we combine 
Eqs. ([2]) and ^ to cancel the space charge term V • E, 
and only then impose the charge neutrality condition 
5p = 6n, where Sn and Sp are the deviations of the 
corresponding densities from equilibrium and we assume 
(|V(5n|/n < 1). The result is 



dt{Sn) = -/iaEext • V(5n + DaV^Sn , 



(4) 



where the e-h grating mobility and diffusion constant are 
defined by 



Ma 



[<7eidn +dp)ah- ahjdn +dp)(7e 
e{(Je + (Th) 

{ae+<Jh) 



(5) 
(6) 



Here Dh = Dhh + Due and = Deh + Dee- By taking 
into account the electron-hole Coulomb drag, the con- 



(PePh - 
^Ph 



PehPh^ + PehPeJ^) ^ 



Peh% 

-Peh 



Pe 



~Peh 



(7) 



where the resistivities of electrons and holes are defined 
as Pe — l/(n/iee) and ph = l/{pphe), respectively (these 
resistivities can be straightforwardly measured in homo- 
geneous d.c. transport experiments in which the carriers 
are either electrons or holes), peh stands for the cross- 
resistivity due to electron-hole Coulomb drag. By sub- 
stituting the conductivity matrix elements into Eqs. ^ 
and ([6]), we obtain 



Pa = 



-^ie{n-p)/n 



1 - 



Pe 



1 + {pe/Ph){Pe/Peh) + {p/n){pe/Ph) {_ Peh 
1 + epehipPe + nph)-enp{ne + p,h){dn + dp)Peh) 



1+Pe/ph + epehPh{n-pY /n 

(8) 

DePtlh + DhU^le + epehPePh{n - p){Dhn - Dep) 



pph + npe + epehPePh{n - pf 



(9) 



These two equations are the main results of this paper. 
For orientation we now discuss some limiting cases, tak- 
ing, for definiteness, the n-type case, i.e., n > p: 
i) Weak Coulomb drag limit {aeh , Peh ~ 0): 



Pa 
Da 



n — p 



n/ph + p/pe ' 

DePPh + DhUpe 

pph + npe 



(10) 

(11) 



These results agree with the well-known expressions for 
the e-h grating mobility and diffusion constant in the ab- 
sence of Coulomb drag Is^] . The e-h grating mobility 
depends on the difference of electrons and hole densities. 
In the n-type case, pa is positive: therefore, the grating 
moves in the direction of the electric field. 
ii) Weak pumping regime {Sn <^ n, ph '3> Pe)'- 



Pa 
Da 



D 



_A P£ 

Peh J 1 + {pe/ph){Pe/Peh)' 

ksT PePe/Peh 



1 + {Pel Ph){Pel Peh) ' 



(12) 
(13) 



where the minority carriers are assumed to be non- 
degenerate. Similar to Eq. (3) of Ref.[28l the expres- 
sion of the mobility here predicts a sign reversal of pa 
at Pe = Peh- Specifically, pa becomes negative when the 
Coulomb drag resistivity is larger than the ordinary resis- 
tivity of the electrons. One notices that the ratio Da/ Pa 
in this case is only dependent on Pe/ Peh and temperature, 
which suggests a neat way to extract the drag resistivity 
in experiment 27 [. 
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Hi) Strong pumping case, i.e., 5n ^ n and, conse- 
quently, p ^ n. In this case, we find that the expression 
in square bracket in Eq. (|5]) is larger than 1 (the derivative 
of Coulomb drag with respect to the density is negative), 
which suggests that Ha can be positive even for pe < Peh- 
If the carrier density due to the injection is much larger 
than the one from doping, then one has n w p, resulting 
in Ha ~ 0. 

In order to test the predictive power of Eqs. ([8]) and 
^ we need a reasonable model for the Coulomb drag 
resistivity peh and the homogeneous drag-free mobilities 
= e/K*E.(^f)-'] andM/. = e/K E.(t')-'], where 
m* and are the effective masses of electrons and holes. 
In two-dimensional polar materials the dominant scatter- 
ing mechanisms are [s^ (i) polar-optical phonon scatter- 
ing ^ 



Aire 



1 e'^Try/2m*oJo/fiNLo ( 1 



2h 



1 

Ko 



, (14) 



where the temperature dependence comes from the 
phonon number N^n, (ii) acoustic phonon deformation 
potential scattering [SJI 



'D 



■E?kBTm* 3 



(15) 



and (iii) ionized impurity scattering, which, in the ran- 
dom phase approximation, has the form 

X (f^ + ^2ekil~ cose)/ Ef) , (16) 
where the screening wave vector is kd = 2toK.o I"^* (1 ~ 

Fermi energy (wave vector), kq and Kqo in these equa- 
tions stand for the dielectric function in the static and 
high frequency limits, respectively, eg is the permittiv- 
ity in vacuum and huQ is the optical-phonon energy. In 
Eq. , D and Vgi represent the crystal density and lon- 
gitudinal sound velocity, respectively, a is the effective 
well-width. S is the deformation potential constant. Fi- 
nally, the resistivity due to electron-hole Coulomb drag 
is obtained from UM 



Peh 



dq- 



e^nphsT (27r)- jq 
X [{l + KD/qfsiiih^{huj/2kBT)] 



dujliaQoelrnQoh 

-1 



(17) 



where Qoe/h is the temperature-dependent Lindhard 
function for electrons or holes times the Fourier trans- 
form of the Coulomb interaction (see Ref. Isil ) . By taking 
the electron doping density uq = 1.9 x 10^^ cm~^ and fit- 
ting the low temperature electron mobility measured in 
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FIG. 1: (Color online) (a) pe (red solid curve) and p^h (blue 
dashed one) as functions of temperature with rn = 0.005no. 
The green dotted curve represents pe with rii = 0.2no. Inset: 
zoom in on the low temperature region T < 3.5 K. (b) Grating 
mobility as function of temperature. The blue squares are 
experimental data from Ref.[23. The apparent breaks in the 
curves as they cross zero are artifacts of the logarithmic scale. 



the experiment by Yang et al. \2'T\ in 9 nm GaAs quan- 
tum wells, /ie(5K) — 5.5 x 10^ cm^/V, we determine the 
effective impurity density to be Ui = 0.005rto. Unfortu- 
nately, the same model, when applied to holes, leads to 
a hole mobility about five times larger than the one im- 
plied by the experimentally measured values of /i^ and 
Da- This discrepancy may be attributed to the band 
structure or additional scattering mechanisms for holes. 
We have found that rescaling by a factor 0.2 at all 
temperatures leads to results in excellent agreement with 
experiment. The following calculations are based on this 
rescaling. 

In Fig.[lja), we plot the electron resistivity pe as a 
function of temperature, as well as the electron-hole drag 
resistivity peh in the weak pumping regime. It is seen 
that the electron-hole drag resistivity is smaller than the 
electron resistivity both at high and low temperature. 
However, at intermediate temperature, one has peh > 
Pe- The two curves intersect at 1.5 K and 385 K. For 
comparison, we also plot the electron resistivity (green 
dotted curve) for a more strongly disordered system, in 
which Hi = 0.2no. In this case, the regime in which 
Peh > Pe is significantly narrower, due to the increase of 
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FIG. 2: (Color online) Grating mobility as a function of 
pumping intensity (the latter quantified by the ratio p/no) 
at 20 and 200 K. 

electron resistivity, but still clearly visible (the electron- 
hole Coulomb drag is essentially independent of impurity 
concentration) . 

Fig-ffl^b) shows the spin-grating mobility fia calculated 
from Eq. ([8]) in the infinitesimal pumping regime. Two 
changes in sign of fia vs temperature are visible: they 
correspond to the two crossings between peh and pe- In 
this figure, we also plot the experimentally determined pa 
(blue squares) . We find that the theoretical results are in 
excellent agreement with the experiment in the temper- 
ature range 20-150 K, in which range the e-h drag domi- 
nates and the e-h grating mobility is negative. However, 
the experimental data do not show any sign reversal. On 
the high temperature end, it simply appears that the 
experiment did not reach sufficiently high temperature. 
A high-temperature cross-over was, in fact, observed in 
Ref. |28|. At low temperature, the reason for the discrep- 
ancy may lie in an undesired laser-induced heating of the 
electrons. 

In Fig. [21 we plot the e-h grating mobility as a func- 
tion of pumping intensity at 20 K and 200 K. Neglect- 
ing the equilibrium hole density, we have p = 6n and 
n — no + p. Then, the pumping intensity can be quan- 
tified by the ratio p/riQ. We find that the mobility at 
first increases with increasing pumping intensity and be- 
comes positive at large excitation levels. After reaching 
the maximum value, the mobility begins to decrease upon 
further increase of the pumping intensity, and eventually 
approaches zero. This behavior is consistent with our 
discussion above. 

Another way to change the sign of the e-h grating mo- 
bility is by tuning the equilibrium electron density by 
means of electrostatic gating. In Fig. [31 we show p,a vs 
Ho for different temperatures and pumping intensities. 
In all cases, the grating mobility at first increases with 
decreasing background density and reaches a maximal 
value. This can be understood by the quick increase of 
the electron resistivity, since Pe/peh increases. As the 
background electron density further decreases to a negli- 
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FIG. 3: (Color online) Grating mobility as a function of equi- 
librium electron density. 
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FIG. 4: (Color online) Ambipolar diffusion constant (blue 
dashed curve) and hole diffusion constant (red solid curve) 
as function of temperature for nt — no in the weak pumping 
case. Notice that the two curves are essentially identical. The 
filled squares represent the experimental data for the ambipo- 
lar diffusion constant reported by Yang et al. 271*] . The green 
dotted curve shows the ratio between the ambipolar diffusion 
constants with and without Coulomb drag, while the pink 
chain curve shows the corresponding ratio for the e-h grating 
mobility. 



gible value, we fall back into the "strong pumping limit" 
{n ~ p) and the mobility decreases towards zero. 

As a last point, we observe that the e-h grating diffu- 
sion constant remains essentially identical to the hole dif- 
fusion constant Dh (see Fig.Hj), both being substantially 
reduced by Coulomb drag caused by interaction with the 
background electrons. We notice that, as expected, the 
effect of Coulomb drag on the grating diffusion constant 
is much less dramatic than its effect on the mobility. The 
theoretical results show good agreement with the exper- 
iment. 

In summary, we have developed a drift-diffusion theory 
for the calculation of the mobility and diffusion grating 
of an electron- hole grating in a single electronic layer, 
taking fully into account the effect of Coulomb drag be- 
tween electrons and holes. Our formulas show that the 
mobility of an electron-hole grating is a sensitive probe 
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of electron-hole drag. Further, due to its dependence 
of electron-hole drag, the mobility of the grating can be 
driven through changes of sign by changing temperature, 
excitation power, or background electron density. 
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